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ABSTRACT 



Aims. To study the Herbig-Haro flows in the nearby dark globule B335. To find new HH objects and H2 knots, make a proper motion 
map of the flow activity and investigate physical properties through shock models. 

Methods. We have observed optical (Ho- and [SII]) and near-IR (2. 12;um Ht) deep fields and taken optical spectra using the 
2.56m Nordic Optical Telescope, as well as a near-UV deep field (if band) using the 3.58 m NTT. In addition we present new 
SPITZER/IRAC (3.5-8.0 yum) and MIPS (24yum) observations. We use previous Ha and 2.12yum H2 observations taken 15 and 9 
years earlier to make proper motion maps. We then investigate the shock physics by matching our spectra with planar shock models. 

Results. We discover five new HH objects (HH 1 19 D-H) in the eastern and one (HH 1 19 1) in the western lobe of the outflow. From 
proper motions we find an optically bright, roughly E~W oriented group with high space velocities (200-280 km s"') and a near-IR 
bright, slower group (15-75 km s"') moving to the ESE. We also find a system of at least 15 H2 knots in the western lobe. This (WNW) 
counterflow suggests the possibility of a binary outflow source, giving rise to two outflow axes with slightly different orientations. We 
find that the E-W flow is symmetrical with evidence for two outbursts. We make the first detection of [01] AA 6300/63 in HH 1 19 B 
and HyS in HH 1 19 A and B and find their extinctions to be Ay ~ 1.4 and 4.4, respectively. HH 1 19 A is found to expand much faster 
than expected from linear expansion with distance from the outflow source. Using planar shock models we find shock velocities of 
~60km s"' (A) and ~35 km s"' (B and C). This agrees with A being of higher excitation than B and C. In our U image we detect three 
of the HH objects and propose that the emission arise from the [Oil] /i3728 line and the blue continuum. New SPITZER/IRAC and 
MIPS observations show most of the HH objects at 4.5 yum and a E-W elongated hour-glass shaped structure at the outflow source. 
Even at 24 fim it is not clear whether most of the light is direct or reflected. 
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1. Introduction 

The nearby prototypical globule B335 is located at about 250 pc, 
as estimated from star counts (Tomita et al. I1979i ). It is very 
isolated and contains a deeply embedded, low-luminosity, low- 
temperature outflow source (Frerking & Langer [T982| l. The ex- 
tinction towards this central outflow source is extremely large 
(Ay > 320 mag, e.g. Chandler et al. 1990). This explains why 
the central object has only been seen at far-IR (Keene et al. 
|1983> and longer wavelengths. It has been observed in sub-mm 
(Chandler et al. 1 19901 1 and at 3.6 cm using the VLA with increas- 
ing resolutions, Anglada et al. 1992 (config D, ±9"), Avila et 
al. [2001 (config C, + 3") and Reipurth et al. ' 2002] (config A, 
sub-arcsec resolution). These observations have shown that the 
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source is variable over timescales of years (Avila et al. 1200 11 1. 
Reipurth et al. I2002l found the outflow source to be clearly elon- 
gated at 3.6 cm in the east-west direction, along the outflow axis, 
indicating a compact thermal radio jet. The combined observa- 
tions makes the B335 outflow source one of the clearest exam- 
ples of a so-called Class source (Andre[T993), the earliest evo- 
lutionary phase of star formation, in which a deeply embedded 
protostar and a bipolar outflow exists. 

The central ~1 x 1 ' has previously been observed in deep H 
and K band images using both Keck I (Hodapp 1998 ) and the 
HST (Harvey et al. 120011 ). Observational evidence for current 
protostellar collapse also exists, based on millimetre wave inter- 
ferometry (Harvey et al. 120031 ). 

B335 is ideal for a case study of isolated star formation, as 
there is absolutely no source confusion caused by OB stars or 
outflows from other young stars, unlike most star forming re- 
gions. Its geometry is similar to that of the well-known HH 30 
(disk and outflow), being seen almost exactly from the side it 
also contains a bipolar molecular outflow accompanied by sev- 
eral Herbig-Haro (HH) objects lying in a line with roughly con- 
stant declination through the deeply embedded outflow source. 
The axis of this HH flow coincides with that of the molecular 
outflow. However, as can been seen in an optical deep field of 
B335 (e.g. the one in this paper) the extinction drops off' rapidly 
about one arcminute from the centre, making it possible to do 
optical spectroscopy of several of the HH objects. 
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The three HH objects known prior to this paper are called 
HH 119B (discovered as a HH candidate by Vrba et al. [19861 
and confirmed by Re ipurth et al. 1992) and HH119A and C 
(Reipurth et al. I1992I I. Hodapp et al. (1998 ) discovered several 
H2 features (called HH 1 19 IR 1-5) using the 2. 12yL(m line of H2. 
In this paper we contribute with a detailed investigation of the 
B335 HH flows using mid-IR Spitzer data, deep near-UV, optical 
and near-IR imaging and optical spectra, revealing six new HH 
objects (two with previously seen near-IR counterparts). We cal- 
culate the proper motion of HH objects and H2 features and for 
HH 1 19 A-C we use planar shock models to calculate shock ve- 
locities, shock compression, preshock densities and extinctions. 

2. Observations and reductions 

2.1. Ground based 

All our ground based observations, except the near-UV imag- 
ing, were carried out at the 2.56 m Nordic Optical Telescope 
(NOT) located at 2382 m above sea level on the island of 
La Palma, Canary Islands. We have used the two cameras 
ALFOSC (Optical spectroscopy and imaging) and NOTCam 
(Near-IR imaging). The near-UV deep field was obtained at the 
3.58 m New Technology Telescope (NTT) located at La Silla 
Observatory, Chile. For details, see the subsections below. 

2.LL Optical spectra 

The ALFOSC (Andalucia Faint Object Spectrograph and 
Camera) long-slit spectra along the optical Herbig-Haro flow of 
HH 119 in B335 were obtained 16 May 2003 using grism #7 
(3850-6850 A) with a slit width of 2'.'5 (R = 260). This slit 
width was suitable since the centres of HH 119 A-C have the 
same declinations within 1'.'3 of each other, however, HHl 19 A 
is extended but with a condensation to the south (HHl 19 A2; 
Reipurth et al. ,1992) that was located inside the slit. With a slit 
length of 6^5 it was possible to observe all three objects in one 
slit for a total exposure time of 2400 s. Even though the Moon 
was full, this did not influence the spectra noticeably since a to- 
tal lunar eclipse was in fact occurring at the time of the obser- 
vations. ALFOSC has a 2048 x 2048 array with a pixel size of 
0'.'188 and the seeing was about I'.'O during the observations. 
There were only three known HH objects (A-C) at the time of 
the observations, but a fourth HH object (D), east of object C, 
was also within the slit due to the alignment of the slit along 
constant declination. 

2.L2. Near-IR deep field (2.12yum H2) 

We have made deep near-IR observations of B335 using the 
2.12jt/m V = 1 - S(l) line of H2. In order to use as much 
telescope time as possible for the line observations, the corre- 
sponding continuum observations were made using a Ks filter 
(2.14/vm central wavelength). The observations were made dur- 
ing 4 photometric nights, 4-8 July 2004, under seeing conditions 
that varied from poor to excellent. We used nb exposure times 
of 119 s and Ks (continuum) exposures of 32 s to avoid satura- 
tion in most bright stars. These exposure times had been tested 
on a previous run and were shown to keep the near-IR camera, 
NOTCam, within its linear regime throughout most of the im- 
age, without any large transient effects from bright sources that 
could otherwise have made a strong pattern in the mosaic from 
the dithering involved. After careful selection of all exposures, 
based on seeing, the final mosaic is both deep and sharp, with 



a mean seeing of 0'.'68 and total central mosaic exposure times 
of 23 205 s (6.45 hours) and 2 880 s (48 min) for the l.ll^im line 
and Ks band respectively. 

NOTCam is an HAWAII 1024 x 1024 xl8.5yum pixel 
HgCdTe array with a field-of-view of 4!0x4'0 (pixel size of 
about 0'.'23). Since the observations are used for proper mo- 
tion calculations we applied a distortion model for NOTCam 
(Galfalk 2005 ) using only the very sharpest exposures by map- 
ping a relation between the on-CCD pixel positions of all known 
2MASS stars in the field to their known coordinates. Using this 
model we have successfully removed most of the previously seen 
distortion, and even at the mosaic edges the PSF of stars is cir- 
cular. 

A median sky is subtracted from the target (with equal ex- 
posure time) and in these observations, where we do not have 
any really extended objects, we use small-step dithering be- 
tween each exposure and calculate a median sky from the on- 
frames themselves. This effectively more than doubles the on- 
target telescope time when compared to beam switching. The 
flatfielding is also differential (skyflats observed with some time 
delay and subtracted). Besides the usual reduction steps of near- 
IR imaging, using in-house routines (written in IDL), we also 
wrote subroutines to find and remove bad pixels, remove im- 
age distortion, sub-pixel shifting and adding overlapping images 
to a mosaic taking into account pixel weighting of individual 
pixel exposure times. We also made a routine to remove any 
dark stripes that results from lowered sensitivity after a bright 
source has been read out of the CCD - this could go unnoticed 
for normal imaging but when a deep field is made with a lot of 
overlaps this can introduce a complicated pattern (especially af- 
ter distortion correction) that has to be removed in order to be 
able to keep a high contrast throughout the mosaic. 

The two mosaics were then sub-pixel shifted in order to co- 
incide, scaled in brightness relative to each other and a con- 
stant added in order to emphasize H2 line emission by having 
the same background levels in both filter mosaics and setting a 
suitable contrast. The final image is a colour composite made 
by RGB colour coding the line image into the red and green 
channels and the (continuum) image into the blue channel, 
resulting in a deep image (Fig. [1]) where H2 emission features 
are clearly yellow while stars, galaxies and reflection nebulae 
are made white/blue. 

2.L3. Western near-IR field (2.12yum H2) 

In May 2007 we made additional S(l) observations, also using 
NOTCam, centred on HH 1 19 A which had not been imaged be- 
fore at near-IR wavelengths. This was also an attempt to find 
new H2 knots in the western outflow. The observations were car- 
ried out on May 28 with total on-target exposure times of 90 
minutes in S(l) and 19.7 minutes in K^. The reductions were 
made in a similar way as desribed in the S(l) deep field sec- 
tion (cf Sect. l2T2] i. The seeing was exceflent (0.31-a.'45 in Ks) 
throughout the night. Part of the resulting image can be seen in 
Fig. |2] which shows three near-IR counterparts of HH 1 19 A and 
a new system of at least 15 S(l) knots to the west. 

2.1 .4. Optical deep field (Ha and [Sll]) 

As a follow up to the near-IR nb-imaging we have also made 
deep optical nb-observations in Ho' and [SII] 6717/31. In addi- 
tion we have made corresponding R band observations to sam- 
ple the continuum for both nb filters. The observations were 
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Fig. 1. Near-IR deep field of B335 in the 2. 12yum S(l) line of H2 (yellow) and Ks (blue) representing the continuum. All previously 
known H2 features are confirmed to belong to the bipolar outflow and their proper motions are shown by the arrows. The red arrow 
shows the high proper motion of a star (asterisk). The position of the outflow source as seen by VLA in the 3.6 cm continuum is 
marked with a plus sign. The proper motions are divided into two groups, HH 119 IR 1 and 5 moving along constant declination 
(similar to the optical HH flow) and IR 2-4 moving more to the south east. However, the groups emanate from the same outflow 
source, as suggested by their proper motions and the cone-like feature at IR 2, with IR 2-4 and IR 5 being on opposite sides of a 
cone shaped cavity created by the outflow. The field of view is 3'89 x 4'13. 



obtained under good seeing conditions (average seeing 0'.'83) 
during 3 photometric nights, 3-6 Aug. 2005 with ALFOSC 
(Andalucia Faint Object Spectrograph and Camera). This instru- 
ment has a 2048x2048 CCD and at a PFOV of 0'.' 188/pixel it has 
a FOV of about 6'4 x 6'4. For both nb filters we used an exposure 



time of 600 s and for the R band 60 s was used in order to avoid 
too much saturation of bright stars. The resulting three mosaics 
have total central exposure times of 24000 s (6.67 hours) Ha, 
22200 s (6.17 hours) [SII] and 2 340 s (39min) for the R band 
(continuum). 
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Fig. 2. Near-IR image of western B335 in the 2.12 fim S(l) line of H2 (yellow) and K5 (blue) representing the continuum. A near-IR 
counterpart of HH 119 A (which we call IR7) is found as well as an extended new system of S(l) features further to the west. The 
HH 1 19 A counterpart consists of three parts, two fairly bright S(l) features (IR7a andb) reminiscent of the optical shock but clearly 
broken up and a fainter part (7c) coinciding with a region of extended [SII] emission in the optical. We find 15 S(l) knots (IR 8a-o) 
in the extended S(l) feature to the west, one (IR8n) which is also seen in Ha and [SII] (Fig.|4]i and classified as an Herbig-Haro 
object (HH 1191). The field of view is 2.'30x 1.'58. 



In addition to the usual reduction steps of optical imaging 
(bias removal, flatfielding, cosmic ray removal) using our own 
set of in-house IDL routines, we also made a camera distor- 
tion model for ALFOSC using our sharpest exposure and the 
positions of 163 USNO-A2 calibration stars. After distortion re- 
moval, all exposures in each filter were sub-pixel shifted and 
combined, resulting in three mosaics. The final colour compos- 
ite mosaic was made with the RGB colour coding Ha (red chan- 
nel), [SII] (green channel) and R band (blue channel). Efforts 
were made to achieve a high contrast in the final mosaic, with 
the same background level in all channels outside of nebulous 
regions and a colour scaling that makes most stars white, nebu- 
lae blue. Ha and [SII] emission features clearly red and green, 
respectively. 

The resulting colour mosiac is shown in Figures[3]and|4] All 
the detected Herbig-Haro objects, except HH 1191 (FigHjl, are 
shown in more detail in Fig.|5] 

2.1.5. Near-UV deep field 

Using EMMI on the 3.58 m New Technology Telescope (ESO) 
during four nights, 28 Jun-01 Jul 2006, we have obtained 
37 x720s exposures of B335 in the U band, resulting in a to- 



tal exposure time of 26640 s (7.4 hours). The filter, U602, is a 
Bessel U band filter centred on 3542 A(covering roughly 3100- 
4000 A) and its transmission curve is shown in Fig.|6l The usual 
reduction steps of optical reductions were applied, including 
bias removal, flat fielding and shift-and-add to correct for track- 
ing errors and small-step dithering between images. The sigma- 
clipping technique was used to remove all cosmic rays after 
making the sky mean value equal in all exposures. 

The EMMI-Blue/NTT 1024 x 1024 chip gives a pixel and 
image field of view of 0'.'362 and 6'.2 x 6'2, respectively. The 
average airmass in these observations is 1.35 and the average 
seeing 2'.' 03. 

2.2. Spitzer Space Telescope 

Spitzer is an IR space telescope that carries a 85 centimeter 
cryogenic telescope and three cryogenic ally cooled science in- 
struments, one of these is the Infrared Array Camera (IRAC) 
that provides simultaneous 5!2x5!2 images in four channels, 
centred at 3.6, 4.5, 5.8 and 8.0 yum. Each channel is equipped 
with a 256 x 256 pixel detector array with a pixel size of 
about 1'.'2 x r.'2. Another instrument is the Multiband Imaging 
Photometer for Spitzer (MIPS) that contains three separate de- 
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Fig. 3. Optical deep field of B335 using a composite of Ha (red), S[II] (green) and the R band (blue). The first thing we note is that 
the molecular cloud surrounding the globule shines unevenly across the image, with both nebulous patches and hole-like features. 
This is probably caused by scattered starlight (in the R band). HH 119 A-E are located in an almost perfect line along constant 
declination, with proper motions showing them to be part of the same bipolar HH flow. Note that the proper motion arrow for object 
F is displaced to fit inside the Figure. The Spitzer contours show the outflow source region at S.Oyum (red, solid curves) and 24/im 
(yellow, dotted curves). The VLA outflow source position is marked with a white plus sign. At 8.0/im there is a clear bipolar-like 
feature caused by reflected light whereas the 24 yum observations have a single peak with a location close to the VLA source. The 
red cross marks the origin HH 119 A would have had based on its current expansion alone (FigfTSli. suggesting that this bow shock 
is currently expanding much faster than before. The field of view is 5fl8 x 4'.2Q. 



Table 1. Spitzer archive data used 



Prog.ID 


P.I. 


Key 


Type 


Released 


Scheduled 


94 


Lawrence C. 


4926464 


IRAC Mapping 


2005-06-24 


2004-04-21 


53 


Rieke G. 


12022016 


MIPS Scan Map 


2005-11-09 


2004-10-15 



lector arrays, making simultaneous observations possible at 24, 
70, and 160 fim. 

The Spitzer data used in this paper (see Table [U were ob- 
tained from the Spitzer Science Archive using the Leopard soft- 
ware. All data had been reduced to the Post-Basic Calibrated 
Data (pbcd) level. We use all channels of the IRAC observations, 
and some in-house routines to reduce them further and especially 
to make a colour composite image. From the MIPS data, we only 



use the 24 /im observations because of its higher resolution and 
because the other two channels are full of artefacts (probably due 
to saturation). The 24 //m camera has a format of 128x128 pixels 
and a pixel size of 2'.'55. In order to plot Spitzer contours on our 
ground-based observations, although very different resolutions, 
we matched all stars seen in the Spitzer 24 jjm mosaic with both 
our Ks and R band observations and corrected for differences in 
FOV, distortion and image rotation. 
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Fig. 4. The westernmost part of our Ha (red), S [II] (green) and R band (blue) deep field. The number of exposures and therefore S/N 
decreases to the west because of flat field uncertainties and image dithering during the observations. After correcting for background 
variations most of the extended system of S(l) knots is covered, however, only one of these bright S(l) knots is detected in the 
optical. This object, IR8n, is seen in both Ha and [SII] and confirmed as a bona-fide HH object (HH 1191). The field of view is 

1:95 X im 



2.3. Proper motions 

2.3.1 . Near-IR first epoch image 

In order for us to map the IR proper motions of as many H2 
features as possible we needed a previous epoch of 2.12 yum ob- 
servations, having a long enough time span to show the proper 
motions with reasonable accuracy, but at the same time with 
high enough quality to show as many of our recently observed 
H2 objects as possible given that the new epoch is a deep field. 
We would like to thank Klaus-Werner Hodapp for kindly pro- 
viding us with the first epoch observations, in the form of a 
deep (8200 s or about 2.3 hours) mosaic taken with the 2.2 m 
University of Hawaii telescope. These observations were made 
using the Quick Infrared Camera (QUIRC) mounted at the f/10 
focus, yielding a pixel size of O'.'ISS (not that far from our sec- 
ond epoch images, 0'.'23). Since these observations were ob- 
tained between 8-10 July 1995 (and our second epoch obser- 
vations 4-8 July 2004) we have a time span of about 9 years 



between the epochs. For more information about the first epoch 
observations see Hodapp et al. 1 19981 



2.3.2. Optical first epocli image 



For our optical proper motions, we used an even larger time span 
(roughly 15 years) than for the IR case. We would like to thank 
Bo Reipurth for kindly providing us with the original Ha im- 
ages from Reipurth et al. 1 19921 taken in August 1990 with the 
3.58 m New Technology Telescope (NTT). These observations 
are deep enough and have a long enough time span when com- 
bined with our images to yield accurate proper motions. The 
instrument used in these observations, EMMI, provided a field 
size of 7.'5 x 7'5 when used with a Thompson 1024 x 1024 pixel 
chip, yielding a pixel size of 0'.'44 (rougly twice the pixel size of 
our Aug 2005 observations). 
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Fig. 5. All Herbig-Haro objects in our optical deep field (HHl 19 A-H), except object I (See Fig|4]l. The colour composition is the 
same as in Fig|3] HHl 19 A and B belong to the west outflow lobe while C-H are part of the eastern lobe. In the case of object C, 
located very close to a bright star, this star has been removed at the position marked with a red plus sign, using the PSF of another 
bright star in the image. The central part of the removed star has been blocked out because of strong residuals that could not be 
removed. Another star can be seen very close to the north of the removed PSF and the C shock itself is fairly extended in both Ha 
and [SII] and shows no strong condensations. Objects A, B, D and F clearly have Ha bow-shocks leading a [SII] cooling zone, 
meaning that the flow moves into a less dense ambient medium as opposed to shock G which moves into a denser medium. Object 
E has a leading shock dominated by a line in the R band that is not Ha or [SII], probably [OI] /l/l 6300/63 emission. The field of 
view in each panel is 38" x 38", except for HH 119 C where the field size is 24" x 24". 



2.3.3. Calculations and accuracy 

We refer to Section 3 in Galfalk & Olofsson 120071 for details on 
our method for calculating proper motions. In short, we fit all 
stars seen in both epoch images using a two-dimensional ellip- 
tical Gaussian. Stars with measurable proper motions are elimi- 
nated in the solution. After warping the lower resolution image 
to fit the one with higher resolution, the images are shifted for 
each object in sub-pixel (0. 1 pixel) steps until the object overlaps 
itself in both epochs, and a proper motion is found. For each step 
the matching is judged by blinking the two images at different 
rates and by making a colour image using two colour channels 
(red and blue). The accuracy is similar for our optical and near- 
IR proper motions (about 12 km/s), except for HH 1 19 C (close 
to a bright star) and HH 1 19 F (faint in the first epoch image). 

3. Results and discussion 

3.1. New HH objects and H2 knots 

Using our deep optical and near-IR mosaics, in combination with 
the proper motions obtained from these and the additional infor- 



mation from optical spectra, we confirm that two of the previ- 
ously known 2.12yL/m H2 features (IR4 and 5 in Hodapp et al. 
119981 1 have optical counterparts that are bona-fide HH objects 
(named HH 119G and D). Four more HH objects are also dis- 
covered in our optical deep field (we call these HH119E, F, 
H and I). These six new HH objects can be seen in Figures[3] 
|4]and|5] HH 119H was first found as an IR H2 feature in the 
Spitzer/IRAC mosaic (FigfTSb. However, it is also faintly de- 
tected in Ha (Fig|5]l and given its position in the E-W HH flow 
it is classified as an HH object. 

Many of the HH objects have the appearance of bow shocks 
with clearly separated Ha and [SII] emission. A summary of all 
shocks is given in Table|2] and the photometry is presented in 
Table |3] Only three of these objects, HH 119 A-C, were known 
HH objects prior to this survey. We have followed the same nam- 
ing scheme, and named Herbig-Haro objects using letters from 
Dtol. 

There is a rough symmetry in HH119A-F and H, lying 
approximately in a line through the outflow source, all at two 
roughly symmetric projected distances, possibly corresponding 
to two outburst episodes in the central source. The innermost 
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Fig. 6. Transmission curves for the UV observations. The three red curves represent the three major components of the system: at- 
mospheric transmission (dotted), detector QE (dashed) and the Bessel U filter (solid). The solid black curve is the total transmission, 
with all three red curves included. Important emission lines in this wavelength range, seen in UV spectra of HH objects (Hartigan et 
al.[l999) are marked by vertical lines. The horizontal line roughly marks the most contributing part of the blue continuum and the 
corresponding spectral energy distribution of the two-photon transition is shown by the blue curve (long dashed) in percentage of 
its peak value. It is clear that only [Oil] 3728 and the blue continuum are likely contributors to the near-UV emission we see from 
three of the HH objects in B335 (A, B and E). 



HH object in the western lobe, HH 119B (58'.'2 from the out- 
flow source), has two counterparts in the eastern lobe, HH 119 C 
(61'.'2) and D (73'.'3). The outermost western object, HH 119 A 
(142'.'3), has three counterparts in the eastern lobe, HH119E 
(123'.'4), F (126'.'4) and H (139'.'4). HH 119G is not located on 
this line roughly along constant declination, instead this object is 
moving more to the south east, in the same group as IR knots 2, 
3, 4a and 4b (Fig.[TJ. All these IR knots are likely to have optical 
counterparts, however, due to the high extinction in this part of 
the flow only knot 4a is detected in the optical and confirmed as 
an HH object. 

In the western lobe we discover a system of 15 new H2 knots 
and three counterparts to HH 1 19 A (Fig.|2]i. 

3.2. Near-IR shocks 

Our 2.12 jjm H2 deep field is presented in Fig[T]along with shock 
names and arrows illustrating the calculated proper motions. The 
near-IR proper motions and projected shock velocities are also 
given in Table|4] There is an apparent division of the shocks into 
two groups, one (IR 1 and 5) that coincides with the optical HH 
flow along almost constant declination, and another group (IR 2- 
4) with a position angle of about 115° (moving to the ESE). 
The projected velocity vectors suggest that both groups emanate 
from the outflow source in the centre of B335 (marked by a plus 
sign in the deep field image). Looking at the east (blue-shifted) 
lobe, the faint cone-like feature close to the outflow source in 
the near-IR (FiglTJ and the mid-IR Spitzer data (FigfTSTi, com- 
bined with the proper motions suggest that the two groups in the 
eastern lobe are on the opposite sides of a cone-shaped cavity, 
cleared out by the molecular outflow. 



While the first group (IR 1 and 5) is bright in the optical 
lines, only the brightest peak in the second group, IR 4a, is seen 
in our optical Ha and [SII] deep fields. The reason for this 
could be higher extinction but another reason is hinted by the 
slower projected velocities, slow moving (C-type) shocks are 
non-dissociative, instead molecules like H2 are kept intact and is 
a good tracer of such shock environments. Fast moving (J-type) 
shocks instead break up molecules, making such shocks bright 
in atomic lines such as Ho' and /l/l 6717, 6731 [SII]. 

West of the outflow source (in the red-shifted lobe) we find 
at least 16 H2 knots, one in the deep field (IR 1) and an extended 
system of knots in the western field (Fig|2]). Although we lack 
proper motions for this system, they are most likely part of a 
counterflow to the ESE flow. This suggests the possibility of a 
binary outflow source, giving rise to both an E-W and an ESE- 
WNW outflow axis. 

The only western knot for which we have previous obser- 
vations, IR 1, shows a proper motion that points almost exactly 
away from the VLA outflow source. Although very diff'erent in 
appearance, this shock coincides very well with the southern part 
of the optically (Ha and [SII]) strong, nicely defined, bow shock 
HH 119B. With both the near-IR and optical shocks having a 
proper motion position angle of roughly 260° and partly over- 
lapping each other, it is clear that they are both part of the same 
HH flow and originate from the central outflow source. However, 
how they are physically related is not equally clear from posi- 
tions and proper motions alone. A similar separation applies for 
the optical and near-IR images of HH 1 19 D (IR 5). A more de- 
tailed discussion in given in Section 13.81 where it is shown that 
both IR 1 and IR 5 likely trace parts of bow shocks HH 1 19 B and 
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Table 2. Geometry of all Herbig-Haro objects and H2 knots. 



HH 119 " 


RA*" 


Dec ^ 


Size^ 


R pcTimp 


^tnir'tiirp ^ 


IVt^l tl iv- c 


f^ommpnt 

V Wl 1 ill it^l 1 1 








\ > 






(first) 




A(IR7) 


19:36:51.32 


07:34:10.5 


6.1 X 17.2 


Optical 


Ha / [SII] 


Reipurth 1992 


Bow shock 


(Ai) 


19:36:51.32 


07:34:10.5 










Ha peak and background star 


(A2) 


19:36:51.47 


07:34:06.8 










Ha peak (South) 


(A3) 


19:36:51.46 


07:34:15.2 




" 






Ha peak (North) 


B (IR 1) 


19:36:56.98 


07:34:06.0 


8.0 X 13.2 


Optical 


Hq'/[SII] 


Vrba 1986 


Bow shock 


C 


19:37:05.00 


07:34:06.0 


1.7x4.4 


Optical 


Ha + [SII] 


Reipurth 1992 


No separation between [SII] and Ha 


D(IR5) 


19:37:05.82 


07:34:07.2 


3.7 X 7.2 


Optical 


Ha/LSII] 


This paper 


Ha bow shock, two [SII] bullets 


E 


19:37:09.19 


07:34:06.5 


5.6 X 8.6 


Optical 


[SIIlVHo- 


This paper 


Leading feature bright in R band"^ 


F 


19:37:09.39 


07:33:46.0 


4.6 X 5.1 


Optical 


Hq'/[SII] 


This paper 


Bow shock 


G(IR4a) 


19:37:04.75 


07:33:38.3 


2.7 X 2.9 


Optical 


[SII] / Ha 


This paper 




H 


19:37:10.22 


07:33:56.8 


5.4 X 4.4 


Optical 


Ha 


This paper 


More extended at 4.5 yum 


I(IR8n) 


19:36:45.41 


07:34:40.2 


2.3 X 3.4 


Optical 


Ha + [SII] 


This paper 


Part of a WNW flow 


IRl 


19:36:57.22 


07:34:01.7 


3.0x5.4 


Near-IR 




Hodapp 1998 


IR counterpart of HH 1 19 B 


IR2 


19:37:01.50 


07:34:08.5 


2.3 X 2.4 


Near-IR 




Hodapp 1998 




IR3 


19:37:03.89 


07:33:49.6 


11.7 X 3.5 


Near-IR 




Hodapp 1998 




IR4 


19:37:05.22 


07:33:39.3 


21.6 X 6.3 


Near-IR 




Hodapp 1998 


Large and irregular 


IR4a 


19:37:04.66 


07:33:36.6 


4.5 X 1.2 


Near-IR 




This paper 


IR counterpart of HH 1 19 G 


IR4b 


19:37:05.77 


07:33:39.3 


3.7 X 2.1 


Near-IR 




This paper 


Condensation within IR 4 


IRS 


19:37:05.70 


07:34:09.5 


13.8 X 11.5 


Near-IR 




T T 1 1 f\r\ 

Hodapp 1998 


IR counterpart of HH 1 19 D 


IR7a 


19:36:51.39 


07:34:12.7 


3.7 X 8.4 


Near-IR 




This paper 


North IR segment of A 


IR7b 


19:36:51.53 


07:34:05.4 


2.8 X 4.7 


Near-IR 




This paper 


South IR segment of A 


IR7c 


19:36:50.93 


07:34:06.5 


2.6 X 3.8 


Near-IR 




This paper 


Segment in front of bow A 


IR8a-o 


19:36:46 


07:34:26 


58 X 33 


Near-IR 




This paper 


Extended system of 15 H2 knots 


IR8f 


19:36:46.23 


07:34:27.2 


0.9 X 1.6 


Near-IR 




This paper 


One of the brightest knots in IR 8 


IR8i 


19:36:46.02 


07:34:29.1 


0.9 X 0.7 


Near-IR 




This paper 


One of the brightest knots in IR 8 


IR8n 


19:36:45.39 


07:34:41.4 


2.3 X 2.6 


Near-IR 




This paper 


IR counterpart of HH 1191 



" Letters designate Herbig-Haro objects, while IR numbers designate H2 knots. 
Optical positions in Ha (04 Aug 2005) and near-IR positions in 2. 12 urn H2 (06 Jul 2004 ; 29 May 2007 for IR 7 and 8). 
Approximate bounding box sizes, determined visually. 
Leading / trailing flux peak. 

The leading feature is easier seen in R than [SII], most likely due to a strong line other than Ha and [SII]. 



D, respectively, that have lower shock velocities (i.e. the outer 
wings of these bow shocks). 

3.3. Optical shocks 

Our optical deep field (Fig.O shows five HH objects 
(HH119A-E) in an almost exact line along constant declina- 
tion. This is a very clear case of a bipolar HH flow, with proper 
motions (Table. |5]i and bow shock morphologies that would be 
expected from their positions relative to the outflow source. 
The projected space velocities of the optically bright and neatly 
lined up HH objects HH119A-D (200-280kms-i) are much 
faster than the near-IR bright SSE moving groups 3 and 4 (15- 
75kms-i). 

One explanation for this could be if the fast moving, east- 
west flow, contains gas expelled into an already cleared out bipo- 
lar cone by the outflow source. The optically bright shocks, HH 
objects, could then be caused by faster gas moving into slower 
gas (with shock velocities much lower than the space veloci- 
ties measured from proper motions). However, at the edge of the 
cone-shaped cavity the expelled gas could interact more with the 
surrounding globule itself via turbulence, making these near-IR 
shocks have slower space velocities, almost stationary, due to the 
encounter with the probably slow moving gas of the globule. 

Given that several HH objects are now known on each side 
of the outflow source, all having similar space velocities, the 
shocks within each lobe have different dynamical ages and there 



has thus been repeated eruptive episodes where matter from the 
accretion is expelled into the outflow lobes. Objects B and C 
are located symmetrically about the outflow source, moving in 
opposite directions, and were probably expelled in the same out- 
burst. Using their (optical) proper motions and distances from 
the central VLA source we find dynamical ages of about 250 
and 265 years for HH 1 19 B and C, respectively. The proper mo- 
tions and dynamical ages of objects B anc C was also measured 
by Reipurth et al. 119921 who suggested that this outburst may 
have occurred 350 years ago. 

There are four more HH objects seen in the optical deep field, 
HH 119F-I. Objects F and H are located in the eastern lobe of 
the E-W flow, close to the HH 119A-E line, while objects G 
(east) and I (west) seem to belong to the tilted (roughly ESE- 
WNW) flow. HH 1191 is the westernmost object in this flow, 
and an optical counterpart of an H2 knot (IR 8n) in a system of 
at least 15 knots. Since these HH objects (and object D) are new 
detections, their proper motions are unknown, except for object 
F which is faintly seen in the Ha first epoch mosaic. Its appear- 
ance, bow shaped with a clear separation of Ha and [SII], and its 
location relative to the outflow source agrees well with the cal- 
culated proper motion. The other objects are clearly HH objects 
even though they lack proper motions, based on their extended 
nature, separation of Ha and [SII], shapes and locations on the 
outflow axis. As a further confirmation of HH status for these 
objects, HH 1 19 E, G and H are clearly seen in the shock sensi- 
tive Spitzer filter centred at 4.5 fim (Fig.fTSl) along with most of 
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Table 3. Flux summary. 
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HH 119 




F [Ha]'' 


F rsiif 


F rSd)!'' 


F4.5 = 


F [SII] / F [Ha]' 












{'->) 




A(IR7) 


8.50 ± 0.34 


17.34 ± 0.91 


9.94 ± 0.47 


17.12 ±0.38 


93 ± 19 


0.573 ±0.018 


B (IR 1) 


1.45 ± 0.07 


05.92 + 0.34 


7.87 ± 0.44 


2.73 ±0.11 


103 ± 13 


1.330 ±0.030 


C 




03.54 ±0.18 


7.36 ± 0.41 






2.079 ± 0.018 


D(IR5) 




01.40 ±0.08 


2.45 ±0.10 


10.95 ±0.17 


107 ± 19 


1.749 ±0.077 


E 


1.23 ± 0.06 


00.79 ± 0.06 


0.51 ±0.03 


out 


40 ±7 


0.65 ±0.11 


F 




00.78 ± 0.05 


2.16 ±0.12 


out 


67 ±8 


2.78 ±0.10 


G(IR4a) 




0.144 ± 0.007 


0.065 ± 0.004 


6.91 ±0.14 


45 ± 11 


0.45 ±0.15 


H(IR6) 




00.07 ± 0.01 




out 


147 ± 17 




I(IR8n) 


out 


01.64 ±0.06 


0.24 ± 0.03 


5.58 ±0.18 


out 


0.148 ±0.021 


IR2 








1.36 ± 0.03 






IR3 








7.71 ± 0.15 


41 ± 6 




IR4 








21.59 ±0.30 


134 ± 13 




IR4b 








2.40 ± 0.05 


14 ±3 




IR7a 








11.04 ±0.34 






IR7b 








5.43 ±0.17 






IR7c 








0.66 ± 0.03 






IR8f 


out 






5.29 ±0.17 


out 




IR8i 


out 






4.34 ±0.14 


out 





" Integrated U band flux in units of 10"'^ erg s"' cm"^, effective filter width is 408.8 A(from total transmission curve). 

Emergent flux in Ha, in units of 10"'^ erg s"' cm"^. 

Emergent flux in the 6717+6731 [SII] lines, in units of 10"'^ erg s"' cm"^. 

Emergent flux in the H2^ = 1-0 S(l) 2.12 jmi line, in units of 10"'^ erg s"' cm"^. 
° Integrated flux of Spitzer IRAC 2 (4.5 ;um) in units of 10"'^ erg s"' cm"^. HH 1 19 C and IR2 close to bright sources. 
' Mean flux ratio including full object extension. 



Table 4. IR proper motion results. Column 1 gives the H2 feature designation, as given in Figure[T] Columns 2-4 give the angular 
distance each shock has moved between the two epochs (9 years apart). Columns 5-7 are the proper motions and finally Columns 
8-9 the (tangential) velocities and positional angles. Uncertainties are given at the bottom of the Table (except for the PA which is 
given for each row). 



HH119 


K 




A,<„ 








vel 


PA" 


Comment 




(") 


(") 


(") 


(mas/yr) 


(mas/yr) 


(mas/yr) 


(km/s) 


(deg) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


IR 1 (B) 


-1.51 


-0.24 


1.53 


-168 


-27 


170 


201 


261 ±3 




IR2 


+ 1.53 


-0.61 


1.65 


+ 170 


-68 


183 


217 


112 ± 3 




IR3 


+0.17 


-0.09 


0.19 


+ 19 


-10 


21 


25 


119 ± 24 




IR4a (G) 


+0.06 


-0.11 


0.12 


+6 


-12 


14 


16 


153 ±37 




IR4b 


+0.45 


-0.35 


0.57 


+50 


-39 


63 


75 


128 ±9 




IR5 (D) 


+ 1.86 





1.86 


+207 





207 


245 


90 ±2 




Star 


+0.17 


-0.35 


0.39 


+ 19 


-39 


44 




155 ± 8 


Marked in Fig.[T] 




±0.06 


±0.06 


±0.09 


± 7 


± 7 


± 10 


± 12 







Position Angle: North ° , East 90 ° , South = 1 80 ° , West 270 ° . 



the other HH objects. The reason for this is discussed in Section 
[3J] 

HH 1 19 E has a different flux distribution than the other HH 
objects. In Fig. [3] its leading shock is seen more clearly in the 
R band than in [SII] emission (making the leading shock look 
blue in the mosaic even though the R and [SII] position of the 
shock overlaps, followed by the Ha shock). This is very likely 
due to unusually bright [OI] AA 6300/63 emission. The situation 
is similar to HH 119B, which has bright [OI] emission lines as 
seen in recent spectra (Fig|2]l; but with an even higher [OI]/[SII] 
ratio since HH 119B looks red-yellow-green across the shock, 
using the same colour scaling. 



HH119G, seen in both Ha and [SII], coincides very well 
with the brightest peak of the corresponding feature in the near- 
IR mosaic (H2 knot IR 4a) and has a separation between Ha and 
[SII] supporting the near-IR proper motion. 



Objects A and E (F) are located roughly symmetrically on 
opposite sides of the outflow source and may belong to the same 
outburst, prior to the B and C outburst. We estimate the ages of 
HH 1 19 A and F to be about 700 and 590 years, respectively. The 
difference of about 100 years in these estimates may be due to 
the large uncertainty in proper motion for object F (which is only 
faintly detected in the first epoch image). 
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Table 5. Optical proper motion results. Column 1 are the HH designations, as given in Figure[3] Columns 2-4 give the angular 
distance each shock has moved between the two epochs (15 years apart). Columns 5-7 are the proper motions and finally Columns 
8-9 the (tangential) velocities and positional angles. Uncertainties are given at the bottom of the table (except for the PA where this 
is given in each row). 



HH119 


A„ 




A,o, 








vel 


PA" 




(") 


(") 


(") 


(mas/yr) 


(mas/yr) 


(mas/yr) 


(km/s) 


(deg) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


A, 


-3.00 + 0.10 


-0.49 ±0.10 


3.04 ± 0.12 


-201 ± 7 


-33 ± 7 


203 ±9 


241 ± 10 


261 ± 3 


A2 


-2.43 ±0.10 


-0.97 ± 0.10 


2.62 ± 0.13 


-162 ± 7 


-65 ± 7 


175 ± 10 


207 ± 1 1 


248 ± 3 


A3 


-2.34 + 0.10 


+0.53 + 0.10 


2.40 ±0.12 


-156 ± 7 


+36 ± 7 


160 ±9 


190 ± 10 


283 ±3 


B 


-3.40 + 0.10 


-0.91 ± 0.10 


3.52 ±0.13 


-227 ± 7 


-61 ± 7 


235 ±9 


278 ± 1 1 


255 ±3 


C 


+3.42 ± 0.19 


-0.57 ±0.19 


3.47 ± 0.22 


+228 ± 13 


-38 ± 13 


231 ± 15 


274 ± 18 


100 ±4 


F 


+3.23 ± 0.57 


-0.57 ± 0.57 


3.28 ± 0.66 


+216 + 38 


-38 ± 38 


219 ± 44 


260 ± 53 


100 ± 10 



Position Angle: North ° , East 90 ° , South = 1 80 ° , West 270 ° . 



Table 6. Optical spectroscopy. 



Object 


Line identification 


Flux inside slit" 


F" 


Fo'^ 


F(.16717)/F(.16731) 








(10-"ergs-'cm-2) 


(Ha) 


(Ha) 




(cm-3) 


HH119A 


H/3^4861 


47.79 ± 4.0 


0.197 ±0.020 


0.304 ± 0.035 








[Ol] ^6300 


20.04 ± 3.7 


0.083 ± 0.017 


0.087 ± 0.019 








[Ol] A6363 


11.73 ± 3.7 


0.048 ± 0.016 


0.050 ± 0.017 








Hff /16563 


242.1 ± 3.6 


1.000 


1.000 








[Nil] ^6583 


10.45 ± 2.6 


0.043 ±0.012 


0.043 ±0.012 








[SlI] ^^67 17, 6731 


75.93 ± 3.7 


0.314 ±0.021 


0.304 ± 0.021 


1.19(1.04-1.36) 


250 (60-500) 


HH119B 


B/3 A4861 


13.3 ±3.3 


0.056 ± 0.016 


0.231 ±0.087 








[Nl] AA5m, 5201 


14.19 ± 2.8 


0.060 ± 0.013 


0.167 ±0.044 








[Ol] /16300 


121.9 ±4.4 


0.517 ±0.027 


0.611 ±0.033 








[Ol] ^6363 


38.35 ± 4.3 


0.163 ±0.021 


0.185 ±0.025 








Hff /16563 


235.8 ± 3.6 


1.000 


1.000 








[Nil] ^6583 


26.27 ± 4.2 


0.111 ±0.020 


0.110 ±0.020 








[SlI] ^^6717, 6731 


440.4 ± 4.2 


1.867 ±0.046 


1.686 ±0.043 


0.82 (0.76-0.90) 


1260 (890-1590) 


HH119C 


[Nl] ^^5198, 5201 


27.64 ± 4.4 


0.138 ±0.028 


-0.384 








[Ol] A6300 


81.46 ± 4.3 


0.408 ± 0.037 


-0.482 








[Ol] ^6363 


35.07 ± 3.7 


0.176 ±0.026 


-0.200 








Ho- A6563 


199.6 ± 7.0 


1.000 


1.000 








[Nil] ^6583 


19.47 ± 5.9 


0.098 ± 0.034 


-0.097 








[SII]^^6717, 6731 


327.3 ± 5.0 


1.639 ±0.085 


-1.481 


1.29 (1.12-1.40) 


140 (30-360) 


HH119D 


[Ol] A6300 


16.19 ±6.8 


0.53 +0.44-0.28 










Ho- .16563 


30.58 ± 6.7 


1.000 










[SlI] ^^6717, 6731 


104.2 ± 6.8 


3.41 +1.24-0.80 









" Using a slit width of 2" 5 along constant declination and centred on HH 1 19 B. 

^ Observed flux, relative to Ha. 

Dereddened flux, using the extinction from the model results found in Table |7] 
Electron density. Estimate of uncertainty range given in parenthesis. 



3.4. Optical spectra and line-ratios of HH 1 19A-D 

The optical spectra of HH 119 A-D are shown in Figures|7] and 
[8] Even though all spectra cover the spectral range 3850-6850 A, 
each spectrum has been plotted using individual ranges, omitting 
parts where no features are seen. 

In Table |6] we present calibrated and relative (to Ha) line 
fluxes for all detected lines. The dereddened flux (Fo) has been 
calculated using the extinctions we find for HH 1 19 A and B us- 
ing planar shock models and comparing the intrinsic and ob- 
served HafH/3 ratios (cf Sect. 13.51 ). We have also calculated the 
mean electron density, n^, in HH 1 19 A, B and C for the position 
covered by the slit using the [SII] I(/16717)/I(/16731) line ratio. 



For this we use the effective collision strengths for electron im- 
pact excitation given in Ramsbottom et al. 119961 and Einstein 
A-coefficients from Keenan et al. 119931 assuming an electron 
preshock temperature of ~ lO'^K which is a common assump- 
tion for HH shock models (see e.g. Morse et al. J 994 ). The 
uncertainty range in these estimates are given in paranthesis in 
Tablei] 

All objects show Ha in strong emission and except for 
HH119A, the [SII]6717/31 doublet is the strongest line. The 
choice of nb filters for our optical deep fields were in fact based 
on the line strengths in these spectra. HH 119 A, B and C also 
have [Ol] AA 6300/63 emission lines ([OI]6300 faintly seen in 
object D as well but with poor S/N). 
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Fig. 7. Optical spectra of HH 119 A and B (NOT/ALFOSC). Object A is clearly of higher excitation than object B. Hfi is detected 
for the first time in object A and B, as well as strong [Ol] lines in HH 119B that were not seen in spectra taken 13 years earlier 
(suggesting variable [Ol] lines). 



It is also very interesting to note that the [Ol] lines are strong 
in HH119B. This is somewhat unexpected, as they were not 
detected in previous spectra by Reipurth et al l 19921 taken 1 3 
years earlier even though the S/N of those observations should 
have been more than enough for a clear detection. This suggests 
that the [Ol] lines are very variable and that conditions have 
changed a lot in HH 119B during this time period. The spectra 
confirm that object A is an intermediate excitation object while 
HH 119B, C and D are found to have lower excitation. Since 
HH 1 19 A is located further out from the outflow source than ob- 
ject B, the lower excitation of B could be explained by it moving 
in the wake of the bow shock from object A. 

It is common practice to plot contours of e.g. Ha-[SII] to il- 
lustrate regions dominated by Ha (positive) and [SII] (negative), 
respectively. This has the disadvantage of cancelling out in re- 



gions where Ha and [SII] have similar flux. It also lowers the 
signal of both the Ha and [SII] contours while the noise remains 
at the same level. This is the case since they always cancel out 
partially as long as there is some Ha or [SII] flux present. A map 
of F([SII])/F(Ha') is better in that respect but has its own disad- 
vantages, the most noticeable one being that very large values 
(both positive and negative) are found where F(Ha) is close to 
zero (mostly in sky regions). It is also common to plot Hq'-i-[SII] 
contours. This uses the full S/N for contours but hides the rel- 
ative flux information that indicates the abundance of the two 
species. 

A solution to this illustrative dilemma is found by e.g. 
using coloured contours to simultaneously plot both the to- 
tal flux F(Hq')-i-F([SII]) (contours) and their relative strength 
F(Ha)/F([SII]) (colours). This is presented for HH119A and 
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Fig. 8. Optical spectra of HH 1 19 C and D (NOT/ALFOSC). These are the first spectra of HH 1 19 C and D, showing both to be of 
low excitation. HH 1 19 C is close to a bright star, therefore some continuum is seen in its spectrum. 



B in Figures |9] and \TU\ respectively. Note that the contours are 
coloured on a pixel-scale level (0'.'188) meaning that each con- 
tour can have any colours along its curve. By choosing reason- 
able surface brightness limits for the faintest contour included in 
the plots only relevant regions are included (i.e. most of the HH 
object and all stars, while excluding the sky). Technically this is 
done by multiplying a high-resolution binary contour map (zero 
for sky and one for contours) pixel-by-pixel with an equally 
high-resolution F(Hq;)/F([SII]) ratio map that has been coloured 
with red and blue shades to separate Ha (blue) and [SII] (red) 
bright regions. 

From these Figures it is clear that HH 119 A is of relatively 
high excitation compared to HH 1 19 B which is a low excitation 
object. There is also a clear separation in both HH 119 A and 
B, between Ha (leading) and [SII] (trailing) in the direction of 
motion. In both cases the bow shock is of high excitation (low 



[SII]/Ha ratio) and the emission behind the bow shock is of low 
excitation (high [SII]/Ha' ratio). This structure is typical for most 
HH objects. Ha tracing sharp leading edges of the cooling zone 
and [SII] dominating in the chaotic and clumpy postshock re- 
gion. It is indicative of a heavy flow that rams into a less dense 
ambient medium. The situation is the same for all other HH ob- 
jects in B335, except for E and G which have [SII] leading Ha, 
suggesting a flow that rams into a denser ambient medium. See 
section l378] for more discussion of the individual HH objects. 



3.5. Shock modeling and extinction 

In this section we use the planar shock models of Hartigan 
et al. ( 119941 ) with our electron densities and different line ra- 
tios (from Table |6]l for HH119A-C to calculate the average 
shock velocity (< Vs >), [SII] -weighted ionization fraction 
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Fig. 9. Detailed overview of HH 119 A simultaneously showing a Hq;+[SII] surface brightness map (contours) and the excitation 
structure by colouring the contours based on the [SII]/Hq' ratio. Note that the colour coding is totally independent of the contours, 
this means that any given contour can have any variation of colours ([Sll]/Hff) along its curve. Stars are also included in this plot, 
these however have black contours. HH 119 A is centred so that its brightest position is roughly located at zero offset. The faintest 
Ha+[SII] contour shown has a surface brightness of about 8.4 x 10"'^ erg s"' cm"^ arcsec"^. 



(< / >), preshock density (no), compression (C) and intrinsic 
(Ha'/Hj6)o ratio of each shock. Using the observed (Ha/Hfi) ratio 
we also estimate the extinction to HH 119 A and B by applying 
the optical extinction law described in Cardelli et al. ( 119891 ) for 
Rv = 3.1, and calculating the extinction Ay needed to match the 
intrinsic ratio from the model. 

This is an iterative process since we need to know some of 
the model results, the preshock density no (for the shock models) 
and the (Hq'/Hj6)o ratio (to deredden the line ratios) to apply the 
models to our observations. The preshock density can be found 
from the models if we can estimate the compression, ionization 
< I > and shock velocity Vs to apply to our electron density 



since no = ("c / < I >)/C. The intrinsic (Ha/H/3) ratio can be 
assumed to be ~3 as a first approximation (case B recombina- 
tion, but usually higher than this for HH objects). By making 
a few iterations through the shock models and recalculating the 
extinction and line ratios each time we can get a converging so- 
lution. However, for each object we need initial rough estimates 
of the shock parameters. For this we compare our spectra with 
the examples given in Hartigan et al. |1994> 

HH119A is an intermediate excitation object, reminis- 
cent of the HH47 jet. We therefore start by assuming that 
Vs ~ 30 km s^' , < / > ~ 2 % and C ~ 20 (from Vs and the shock 
models). Using our observationally determined n(,~250cm"^ 
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Fig. 10. Detailed overview of HH 119 B, using the same technique as in Figure |9] HH 119B is centred so that its 
brightest position is roughly located at zero offset. The faintest Hq;+[SII] contour shown has a surface brightness of about 
3.1 X 10"'^ erg s"' cm"^ arcsec"^. 



we estimate that no ~ 625 cm . For a (case B) (Ha/HjS)o ~ 3 and 
the observed value 5.076 the extinction becomes Ay - 1 .65 and 
E(B - y) = 0.53. The dereddened line ratios [OI] 6300/Hq', 
[SII] 6717+31/Hq' and [Nil] 6583/[OI] 6300 then gives three 
estimates of Vs , < / >, C and (HcK/HyS)(), one for each line ratio. 
These values are then used for the next iteration. 

HH 119B and C are of low excitation, similar to the HH34 
jet (Hartigan et al. |1994| l. we therefore make the initial estimates 
Vs ~ 30kms-i, < / > ~ 1 % and C ~ 15. 

The fact that we do not detect [OIII] /15007 in any of the 
objects suggests (Morse et al. 1993) that the shock velocities in- 
volved are lower than 90 km s" ' (for no preionization), 80 km s ' 
(equilibrium preionization) or 40 km s"' (fully preionized). This, 



in addition to the excitation types, supports our initial estimates 
of slow shock speeds (V^ ~30kms"'). 

The initial estimates, iterations and results are given in 
Table |7] Since we do not detect Hj8 in the spectrum of HH 1 19 C 
its extinction is assumed to be roughly the same as for object 
B, motivated by the fact that their projected distances from the 
outflow source are about the same. The resulting extinctions are 
found to beAy = 1.36 (1.05-1.67),£(B-y) = 0.44(0.34-0.54) 
for HH 119 A and Ay =4.44 (3.67-5.28), £(B-y) = 1.43(1.18 
- 1 .70) for HH 1 19 B. The assumption that object C has the same 
extinction as object B based on equal distance from the outflow 
source may however be inaccurate, as B335 has been shown to 
have a long tail (Frerking et al. 1 19871 1, a remanence from its for- 
mation history as a cometary globule (Reipurth il983) . However, 
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Table 7. Extinction and shock model results for HH 119 A, B and C. 



Iteration step 


{Ha/HJ3)o 


A/ 


E(B-V) 


no 


C Ionization I 


< / > = 


Shock velocity V^"* 


< Vs > 






(mag) 


(mag) 


(cm-3) 


(%) 


(%) 


(km/s) 


(km/s) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) (7) 


(8) 


(9) 


(10) 


HH 119A 


Initial guess ' 


3.00 


1.65 


0.53 


625 


20 


2.0 




30.0 


Iteration 


3.34 


1.31 


0.42 


25 


40 40.0, 23.0, 13.5 


25.5 ± 14 


74.7, 60.2, 46.4 


60.4 ± 15 


Result 


3.29 


1.36 


0.44 


23 


40 40.0, 30.0, 13.5 


27.8 ± 14 


74.7, 68.7, 45.3 


62.9 ± 16 


HH 119B 


Initial guess ^ 


3.00 


5.59 


1.80 


8400 


15 


1.0 




30.0 


Iteration 


4.33 


4.44 


1.43 


1530 


25 3.1,2.2,4.6 


3.3 ± 1.3 


37.8, 34.6, 39.8 


37.4 ± 2.7 


Result 


4.33 


4.44 


1.43 


1070 


28 3.7, 3.6, 5.2 


4.2 ± 1.0 


37.1, 35.7, 38.0 


36.9 ± 1.5 


HH 119C 


Initial guess ' 


4.33 


(4.44) s 


(1.43)*? 


930 


15 


1.0 




30.0 


Iteration 


4.18 


(4.44) s 


(1.43) 


100 


28 6.1,4.1,6.1,3.3 


4.9+ 1.5 


43.6, 37.7, 38.1,33.5 


38.2 + 4.2 


Result 


4.38 


(4.44) 




83 


23 8.6, 7.6, 7.7, 5.5 


7.3 ± 1.5 


41.9, 37.2, 33.3, 26.9 


34.8 + 6.4 



" Using tiie extinction law of Cardelli et al. |1989l Rv = 3.1 and (Ha/H^)o6i. 
Preshock density, Hq = (n^/ < / >)/C. 
Shock compression. 

Estimates from line ratios [01] 6300/Hq', [SII] 6717+31/Ha, [Nil] 6583/[OI] 6300. For HH 1 19 C also from [NI] 5200/[NII] 6583. 

[SII]-weighted ionization fraction < I >. 
' Initial guesses taken from objects with similar excitation, HH 47 jet (HH 1 19 A) and HH 34 jet (B and C). See Hartigan et al l 19941 
8 Since HyS is not detected from HH 1 19 C, the extinction is assumed to be the same as for HH 1 19 B (see text). 



since we only use line ratios close to each other in wavelength (a 
few hundred A), except for the [NI] 5200/[NII] 6583 ratio, the 
model results for HH 119 C should not be very sensitive to this 
assumption. 

It is interesting to note that shock velocity and ionization 
found for object A is much higher than that for object B. This 
agrees with the suggestion that object B is following in the wake 
of the outermost object A and thus colliding with cloud material 
(or turbulence in the flow itself) that has already been sped up by 
object A. The higher excitation of HH 1 19 A also agrees with a 
higher shock velocity and ionization. HH 1 19 B is found to have 
much higher preshock density no than object A, which might 
be expected given its location much closer to the centre of the 
globule where the density is higher. 

3.6. Near-UV imaging 

In Figure [TT] our U band deep field is presented. HH119A 
is clearly detected, and has a very similar appearance to its 
Ha counterpart. After background subtraction of the scattered 
starlight in the globule we also detect HH 1 19 B and E but these 
are however much fainter than object A in the near-UV. For ob- 
ject B this is expected because of its much higher extinction 
(AAy w 3). None of the other HH objects are detected, but for 
HH 1 19 C and D this could partly be explained by their proxim- 
ity to a bright star combined with the seeing in the near-UV. 

The transmission curve of the U filter is shown in Fig. [Alto- 
gether with the atmospheric transmission for the mean airmass 
of the observations, the detector QE and the total system trans- 
mission (black solid curve). There are probably only two spec- 
tral features that contribute to the U band flux of these HH ob- 
jects. The [Oil] /13728 line, located close to the maximum sys- 
tem transmission, which has also been shown to be strong in 
other HH objects (see e.g. Hartigan et al. II 9991 ) and the blue con- 
tinuum from two-photon emission. The blue continuum spectral 



distribution has been plotted in Fig. |6] using a blue dashed curve 
with a scale corresponding to the percentage of its peak value. 
The total U band flux for the three detected HH objects are given 
in Table H 

3. 7. Spitzer mid-IR imaging 

FigurefTSlshows our Spitzer 3.5, 4.5 and 8.0 fim IRAC composite 
of B335. An hourglass-shaped reflection nebula is clearly seen 
in the centre (at 4.5 and 8.0jum), with an outflow cone to each 
side of the central VLA source (IRAS 19345-H0727), marked by 
a cross in the Figure. Most of the HH objects we see in our opti- 
cal and near-IR deep fields can also be seen in the 4.5 /urn chan- 
nel of Spitzer (green in the colour composition used in Fig.fTsTl. 
Note that, except for the strong artefacts seen close to bright 
stars (to the SSW) in this filter, all these extended objects are 
most likely HH objects. These include all HH objects A-H, ex- 
cept HH 1 19 C, which is too close to a bright star for detection. 
HH 1 19 1 is located outside the Spitzer image. 

It has been shown in several other Spitzer surveys (e.g. 
Noriega- Crespo et al. |2004l Harvey et al. |2006l and Galfalk & 
Olofsson 120071 ) that the IRAC channel centred at 4.5 fj.m is very 
efficient in detecting bona-fide Herbig-Haro objects. The rea- 
son for this is partly that the spectral response function is high- 
est in this channel, but there are at least two more contributing 
factors. Between approximately 4-5 fim there are many vibra- 
tional and rotational H2 emission lines, these have been mod- 
eled by Smith & Rosen (2005 ) for afl IRAC bands using three- 
dimensional hydrodynamic simulations of molecular jets. The 
strongest integrated H2 emission is predicted to arise from band 
2 because of rotational transitions. For the typical conditions of 
low-mass outflows, pure-rotational transitions like S(ll)-S(4) 
(4.18-8.02jum) can actually be much brighter than the standard 
2.12/im H2 line (Kaufman & Neufeld |1996l l. Channel 2 is also 
the most "PAH-free" band of IRAC, greatly enhancing its use- 
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Fig. 11. Near-UV deep field. A composite of Ks (red), R (green) and U (blue) mosaics (left panel). The cloud extinction is effectively 
shown using this filter combination, spanning from the near-UV to the near-IR. After removing the background from scattered 
starlight in the U mosaic, three of the HH objects are detected (HH 119 A, B and E). In the right panel we show several filter 
combinations using background subtracted mosaics. It is clear that the appearance of object A is very similar in the U and Ho; 
filters. The field of view is 4'.25 x 4f 13 in the left panel and 1 16" x 39" in each of the right panels. 



fulness as a HH tracer, as opposed to the 5.8 and 8.0 //m channels 
which have very high background contrasts caused by extended 
Polycyclic Aromatic Hydrocarbon (PAH) emission, hiding the 
shock-excited H2 features of the HH flows. 

A closer inspection of our IRAC composite actually reveals 
one more extended object with strong 4.5 yum flux, this object 
is located to the east of HH 1 19 E and F, at a declination that is 
consistent with it being part of the HH flow. It is located just out- 
side the eastern border of our 2.12 fim mosaic but a faint object 
is hinted at this position in our Ha mosaic, although very close 
to the border where the S/N is degraded due to the few number 
of overlapping exposures used in that part of the mosaic. Given 
its location and brightness in the Spitzer image, together with its 
Ha detection (although faint), we classify this as an HH object 
(HH119H). 

We identify a very red object to the north in Fig. [15] as 
HH119VLA2 (Avila et al. 1MB, also known as Barn 335 3 
(Anglada et al. |1992| l. The nature of this source is unknown. It 
is located outside the northern border of our 2.12/im H2 S(l) 
mosaic and no Ha emission is seen near this object. 

3.8. Individual Herbig-Haro objects 

Figure|5] shows all the HH objects in detail with individually 
adapted contrast. In this Section details of each HH object will 
be discussed, referring to Fig.|5]in the following if the image of 
a shock is discussed without further reference. 

3.8.1. HH119 Aand IR7 

The image of HH 1 19 A has been deconvolved using maximum 
entropy to bring out more details, made possible by the high S/N 
ratio achieved for this object. The leading bow shock is domi- 



nated by Ha emission (even though there is considerable [SII] 
emission as well) followed by a region dominated by [SII] emis- 
sion. There are three condensations (emission peaks) in the bow 
shock, called A 1 (middle peak), A2 (south) and A3 (north). 

We have measured the proper motions within the bow shock 
using the three Ha peaks. The result is shown in Fig.[T2]with pro- 
jected velocities illustrated by arrows. The image is an extreme 
close-up of the HH 1 19 A part of our deep field, with overplotted 
Ha contours adapted from Fig. 2 of Reipurth et al. 1992 Note 
that even though the northern condensation (A3) may seem to 
move in a slightly different direction (not as far north) than the 
arrow implies in the left panel of Fig. [121 this is just a visual ef- 
fect caused by both Ha and S[II] being shown as red and green, 
respectively, in the colour composite. The proper motions were 
calculated from the locations of the Ha emission knots, as shown 
in the right panel. 

We note that the central peak found in the first epoch image 
coincides with a star behind the shock that is now revealed 15 
years later The Ha contribution from this star in the first epoch 
is unknown, but a wider peak has been used to trace the proper 
motion of the middle part of the shock in order to minimize the 
risk of confusion with this star 

The bow shock is clearly expanding at a much faster rate 
than would be expected by its distance from the outflow source 
assuming a linear expansion, since tracing the arrows back- 
wards leads to a point much further west than the outflow source 
(Figures[3] andfTSI). The expansion rate has thus increased with 
time as the shock gradually moves into less dense parts of the 
globule, approaching free expansion. Object B has a similar ap- 
pearance to object A, albeit on a smaller scale, suggesting the 
amount of expansion between these positions in the flow. The sit- 
uation is reminiscent of the HH34 jet and bow shock (Figures 9 
and 10 of Reipurth et al. 120021 ) where the width of the jet and the 
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Fig. 12. Optical close-up of the HHl 19 A shock in Ha (red), S[II] (green) and R (blue). The overplotted Ha contours from the first 
epoch (15 years earlier, Reipurth et al. 1 19921 1 clearly illustrates the movement of this shock. Proper motions have been calculated 
for three Ha peaks in the shock (Ai, A2 and A3) as shown by the arrows. The shock is clearly expanding much faster than would 
be suggested by its position in relation to the outflow source if a linear expansion was assumed. The field size in the left panel is 
25'.'9x27'.'8. 



size of HH objects increase with distance along the jet from the 
outflow source. 

The spectrum of HHII9A2 (Fig.]?) shows strong Ha and 
weak [SII] emission, HH 119 A is thus a high excitation object 
as opposed to the other HH objects for which we have optical 
spectra. The first detection of H/3 in object A is also seen in this 
spectrum. 

The picture of HH 1 19 A as a coherent bow shock (as seen in 
Ha) is not always true. We detect the S(l) counterpart of object 
A as broken up into three parts (IR 7a-c, see Fig.|2]l. In Figure[T3] 
we plot flux contours in Ha, [Sll] and the S(l) line. It is ev- 
ident that most of the S(l) knots in the northern and southern 
part (IR 7a and b) lie well within the contours of the optical bow 
shock. There is also a third, very faint counterpart in S(l) (lR7c), 
also bow-shaped, which preceeds the optical bow. Comparing 
the [Sll] and S(l) contours (right panel) we see that, not only is 
there extended [Sll] emission in a cooling zone behind the bow 
shock, there is also a faint [Sll] structure preceeding the bow, at 
the position of IR7c. 

3.8.2. HH119 Band IR1 

HH 119B is a nice well-defined bow shock with Ha followed 
by [SII] emission, and an apex with both strong Ha and [SII] 
emission. The shock follows in the wake of HH 1 19 A, which is 
a possible explanation why the spectrum of this shock shows a 
lower excitation (see Fig.|2ll. The spectrum also shows object B 



to have strong [OI] AA 6300/63 emission lines, which was not 
the case in the Aug. 1990 observations of Reipurth et al. ( 19D2| l. 
The conditions in object B have apparently changed during the 
13 years between the two epochs. 

The near-IR H2 S(l) source IR 1 has a similar position and 
proper motion position angle as object B, although their proper 
motions differ by more than 70kms~'. To clarify the situation a 
colour composite using both the optical and near-IR shock ob- 
servations has been made (Fig. [141 right panel) by warping the 
S(l) image to match the pixel size and distortion of the optical 
images (Ha and S[II]). 

These shocks are Ukely physically connected, with the S(l) 
bright shock IR 1 located in the southern part of optically bright 
shock B. While their difference in projected velocities suggest 
that their alignment could be merely coincidental, the most likely 
situation is that they are in fact counterparts of the same object, 
since IR 1 is located far out in the wing of bow shock B where 
the (projected) shock speed is much slower than at the apex. 
This could explain the S(l) brightness if H2 is not disassociated 
in the bow wing. The geometry and expected excitation struc- 
ture (Fig.O thus suggest that IR 1 is the near-IR counterpart of 
HH119B. 

3.8.3. HH119C, Dand IRS 

In the image of object C in Fig|5] a bright star has been removed 
at the position marked by a plus sign in all three filters (Ha, [SII] 
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Fig. 13. HH 1 19 A flux contours in Ha (blue), S[I1] (green) and S(l) (red). The contour levels have been chosen to emphasize the 
different features seen in each transition. It is clear that most S(l) knots in HH 1 19 A lie well within the optical bow shock. In S[ll] 
there are two more features when compared to Ho-, a bright trailing structure and an extended faint feature in front of the bow. This 
feature can also be seen faintly in S(l) as a bow-shaped object leading the main optical bow shock. Stars have been coloured black. 
The field size is 21" x 21" in both panels. 






Fig.l4. Optical and Near-IR composites of HH 119D/IR5 (left) and B/IR 1 (right) using a Ha (blue), [SII] (green) and 2.12//m 
H2 (red) colour coding. In both shocks the optical shock emission is clearly strongest close to the apex of the bow shock. Ho- at the 
front and [SII] in the cooling zone behind it. The near-IR shock emission is displaced to the north and south in HH 119D and B, 
respectively. In the case of HH 1 19 D, the near-IR shock is also bow-shaped, but located further to the north. The near-IR counterpart 
of HH 1 19 B, IR 1, does not have the appearance of a bow, it does however coincide with the far southern wing of the optical bow 
shock, as might be expected since the projected shock velocity (perpendicular to the bow surface) is much lower far out in the wing 
than at the apex. Both fields have a size of 75" x 75". 
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and R) using the PSF of a similarly bright star, lowering the con- 
trast considerably and making it possible to display HH 119 C in 
much more detail. A star, not seen prior to the PSF removal, is 
also revealed very close to the north of the removed bright star 
(at the edge of the black circle). 

This is the first detection of object D in the optical, it had 
only been seen faintly in the near IR prior to this paper (IRS). 
Our near-IR observations show this to be a bow shock at 2.12 fim 
H2 as well. Object D has a complicated excitation structure, as 
suggested by the optical/IR composite image (FigfT4li. A leading 
bow-shock dominated by Ha is clearly seen, followed by two 
well-defined 'sulphur bullets' in the complicated cooling zone 
and a near-IR H2 bow shock, displaced to the north relative to 
the optical bow shock. The situation appears to be similar to that 
of objects B and IR 1, with a S(l) strong part of the bow shock 
far out in the (north) wing where the projected shock speed is 
much lower than at the apex of the bow shock. 

3.8.4. HH119 E, F, H and IR6 

Both HH 1 19 E and F are previously undetected Herbig-Haro ob- 
jects. As can be seen in the image, they both have separated neb- 
ulous features between filters in the direction expected if they 
originate from the central outflow source. HH 119 F is seen as a 
Ha bow shock, followed by a similar S[II] feature. It resembles 
a smaller version of HH 119B moving approximately ESE, as 
suggested by its structure and supported by its proper motion. 
HH 119E is probably moving almost exactly to the east, as it 
is perfectly lined up with HH 119 A-D and the outflow source. 
As seen in the image for object E, and suggested by the optical 
spectra of other HH objects in this flow (B and C), object E must 
have very strong [OI] AA 6300/63 emission (as the leading shock 
is strongest in the R band) followed by a Ha shock. 

There is another, although very faint. Ha feature to the east 
of these objects in the HH flow but mid-way between them in 
declination. This object, HH 119H, is however clearly seen at 
4.5 yum using Spitzer (IR6 in Fig[T5]l. Details about this object 
are given in Table|2] 

3.8.5. HH119G, IR2, 3and4 

In the blue-shifted lobe we find IR objects 2-4 moving as a group 
roughly to the ESE. Except for the brightest S(l) peak, IR4a 
(G), these objects are only seen in the near-IR deep field. They 
all have slower projected velocities than the A-D optical bipolar 
flow, and especially features 3 and 4 have much slower space 
velocities. 

As suggested by the Spitzer image (Fig.flSll these shocks 
could be the result of turbulent interaction of molecular mate- 
rial from the surrounding cloud and the outflow, near the cavity 
walls. This interaction could be what is slowing down the out- 
flow and creating these shocks. The reason why they are only 
seen in the near-IR could be a combination of extinction and the 
amount of H2 relative to atoms and ions (such as H, S+ and O). 
In a slow magnetically supported (C-type) shock, molecules are 
not as dissociated as in fast (J-type) shocks. IR2 is moving much 
faster than objects 3 and 4, in a direction more to the east, within 
the cone that has already been cleared out by the outflow. 

3.9. HH1 19 I and the IR 8 system 

In our 2007 images of western B335 we see a very likely coun- 
terflow (Fig.|2]i to the ESE flow, consisting of a system of at least 



15 H2 knots. Careful image re-reduction of our optical deep field 
close to the mosaic edge, where the signal is lower due to few 
overlaps and the flatielding is very sensitive to vignetting, re- 
vealed an optical (Ha and [SII]) counterpart, HH 1191 (Fig.|4]i, 
to one of these H2 knots. This WNW counterflow (although 
no proper motions are known yet) is, similar to the ESE flow, 
much brighter in S(l) than at Ha or S[II] when compared to 
HH 1 19 A-F in the E-W flow. The detection of this system and 
the similar properties to the ESE flow opens up the possibility 
that the outflow source may in fact be a binary source, with two 
slightly different oriented outflow axes giving rise to these two 
flows. 

4. Summary and conclusions 

In an effort to investigate the Herbig-Haro flows in the nearby 
dark globule B335, a well-studied protostellar collapse candi- 
date with three previously known HH objects, we have taken op- 
tical spectra and observed an optical tri-colour deep field (Ha 
6.7 h, [SII] 6.2 h and R 39min) and a near-IR bi-colour deep 
field (2.12yum H2 S(l) 6.5h and Ks 48min), as well as an ad- 
ditional near-IR field of western B335, using the 2.56 m Nordic 
Optical Telescope and a near-UV deep field (U band 7.4 h) with 
the 3.58 m New Technology Telescope. In addition we present 
new SPITZER/IRAC (3.5, 4.5 and B.OA^m) and MIPS (24/im) 
observations. The following results were obtained: 

1. We discover five new HH objects (HH 119D-H) in the 
blue-shifted (eastern) lobe and one new HH object (HH 1 19 1) in 
the red-shifted (western) lobe of the outflow. Two of these HH 
objects (D and G) have previously observed near-IR counterparts 
(Hodapp |1998l l. In our S(l) mosaic we detect these as large ex- 
tended shocks (object D being a large bow shock). In western 
B335 we also detect a system of at least 15 previously unknown 
H2 knots. 

2. By combining our observations with earlier epochs in the 
optical and near-IR we calculate proper motions of the previ- 
ously known HH objects (and the new HH 119F) and 2.12/im 
H2 S(l) knots using a time span of 15 and 9 years, respectively. 

3. The combined optical and near-IR proper motion maps, 
together with the deep images, show that the opening angle 
of the bipolar flow is already very wide and reveals two sub- 
groups within the flow, a high velocity bipolar flow (space veloc- 
ities 200-280kms-i) in the E-W direction (HH 119 A-F) that 
is bright in optical shock lines, and a much slower flow (15- 
75 km s ' ) in the ESE direction within the blue lobe (HH 119 G 
and IR2^) that is bright in the S(l) line. We also propose a 
counterflow in a roughly WNW direction (red lobe), also bright 
in S(l) and faint in Ha and [SII]. 

4. All proper motion vectors originate from central B335. We 
suggest two possibilities, either the E-W and ESE- WNW shock 
groups belong to the same outflow but are flowing in different 
parts of a cone-shaped cavity cleared out by the outflow, or the 
outflow source is binary in nature resulting in two outflow axes 
with slightly different orientations. 

5. Two of the new HH objects, HH 119E/F, are equidistant 
to HH119A from the outflow source, while HH119C/D are 
roughly equidistant to HH 1 19 B. The bipolar E-W flow is thus 
roughly symmetrical, with two outburst episodes suggested, 255 
and 645 years ago, respectively, from proper motions. The con- 
clusion that the central source has undergone several eruptive 
events has previously been suggested by Reipurth et al. (1992]l 
who also used proper motions to estimate these two outburst 
ages (they estimated 350 and 850 years). 
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6. For HH 1 19 A, the proper motions of three Hor condensa- 
tions show that the shock is rapidly expanding, at a much faster 
rate than would be expected from a linear expansion given its 
distance from the outflow source. 

7. The [SlI] / Ha structure of the shocks indicate that the ESE 
group rams into a denser ambient medium, with the situation 
being the opposite for the E-W flow. 

8. Our optical spectra show HH 119 A to be of higher exci- 
tation than B and C (which are both low excitation objects). We 
detect H/S for the first time in both A and B, making it possi- 
ble to estimate the extinction towards these HH objects. In stark 
contrast to a previous spectrum of HH 119B, we detect strong 
[01] AA 6300/63 emission lines. It is clear that object B must 
have changed a lot during the last decade or so. 

9. We calculate the optical extinction towards HH 1 19 A and 
B to be Ay w 1.4 and 4.4, respectively. 

10. Using planar shock models we have estimated the shock 
properties of HH119A-C. The resulting shock velocities are 
much lower (a factor 5-10) than the space velocities found 
through proper motions. This is however very common for 
Herbig-Haro objects, and expected when gas with high space 
velocity is caught up by slightly faster gas, making the relative 
(shock) velocity lower than the space velocities. This is also the 
case in bow shocks, where the shock velocity is lower in the 
wings than at the apex of the bow since the projected velocity 
perpendicular to the shock is highest at the apex and lowest far 
out in the wings. We find shock velocities of ~60km s ' (A) and 
~35kms ' (B and C). This agrees well with object A being of 
higher excitation, and B following in the wake of this object. 
The [SII]-weighted ionization fraction is also found to be much 
higher in object A. 

1 1 . We detect HH 1 19 A, B and E in the U band after back- 
ground subtraction of scattered starhght throughout the globule. 
This emission is proposed to be a combination of two spectral 
features, the [Oil] A372S line and the blue continuum. 

12. Most of the HH objects and S(l) knots are detected in 
the SPITZER/IRAC observations (channel 2, 4.5 yum). An hour- 
glass shaped feature is seen in the IRAC images, centred on the 
VLA outflow source, with a morphology that agrees with a cone- 
shaped E-W flow. Even at 24 yum it is unclear whether most of 
the light is direct or reflected, a single peak can be seen but it is 
offset by about 2'.' 8 to the SE from the VLA position. 
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Fig. 15. Spitzer IRAC composite of B335 using a 3.5 yum (blue), 4.5 /im (green) and 8.0/im (red) colour coding. The red plus 
sign marks the central VLA source position. There is a clear correlation between the bipolar cone-shaped structure seen at these 
wavelengths (3.5-8.0//m) and the VLA source position. Many of the HH objects and S(l) knots seen in our optical and near-IR 
mosaics are also detected at 4.5 //m (IRAC channel 2) which is clearly the most efficient channel for tracing the flows. 



